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Abstract. The binary LS 5039 is a non-thermal X-ray emitter that presents 
jet-like radio structures, and is also one of the most misterious TeV sources in 
our Galaxy. The presence of an O-type star in LS 5039 implies that the non- 
thermal emitter must be embedded in a strong stellar wind, and the role of the 
latter could be relevant for the understanding of the high-energy behavior of the 
source. In this work, we show that the lack of absorption features in the soft 
X-ray spectrum of LS 5039 can constrain strongly the parameters that describe 
the wind, and ultimately the location of the non-thermal emitter. 



1. Introduction 

The X-ray binary LS 5039 presents non-thermal radio and very high-energy 
emission (e.g. Marti et al. 1998; Aharonian et al. 2005), and the X-rays emitted 
by this object are also of likely non-thermal origin (e.g. Bosch-Ramon et al. 
2007 -B07-; see also Takahashi et al. 2009). The system consists of an O main 
sequence star and a compact object of unknown nature (Casares et al. 2005 
-C05-), and presents jet-like or collimated radio structures (e.g. Paredes et al. 
2000, 2002; Ribo et al. 2008). 

LS 5039 was observed in 2005 by XMM-Newton right after the periastron 
passage (phase ~ 0.04; OBS1), close to the superior conjunction of the compact 
object (SUPC=0.06), and around the apastron passage (phase ~ 0.5; OBS2), 
slightly before the inferior conjunction (INFC=0.72)Q- These observations, pre- 
sented in detail in B07, showed that photo-electric absorption may only take 
place in the interstellar medium, being absorption in the stellar wind negligible, 
and put an upper-limit of 10 21 cm~ 2 on the difference between the absolute hy- 
drogen column densities of OBS1 and OBS2 (AiVn). This upper limit is about 
an order of magnitude smaller than the value predicted adopting a homoge- 
neous spherical stellar wind (B07). Possible explanations put forward in B07 
for such a small photo-electric absorption were a strongly ionized wind, which 
seems unlikely due to the moderate X-ray luminosity of the source of few times 
10 33 erg s , an X-ray emitter far from the compact object at the borders of the 
binary system, and finally an inhomogeneous wind with a large porosity length 
(see, e.g., Owocki & Cohen 2006). Here, we explore the possibility that the 



New ephemeries have been presented by Aragona et al. (2009), although they are similar to 
those of C05 and do not change the conclusions of this work. 
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porosity length is long enough to effectively reduce 7V H at large, and investigate 
what are the implications of this for the stellar wind. 

2. Constraints on the stellar wind properties from soft X-ray data 

The constraints on the Nr inferred by B07 can be used to derive restrictions 
on the properties of the stellar wind of LS 5039. For this, we adopt a simple 
clumpy spherical wind model in which the clumps are assumed to be spherical 
with a constant radius of R c = 10 10 R c io cm. The wind volume filling factor 
is fixed to its value fo = dVo/dVo = 0.01 /o.oi at the periastron distance a = 
ro ~ 1.4 x 10 12 cm from the star. Here, dV ' and dVo are the volume occupied 
by clumps and total one, respectively, for a shell located at ro. At the same 
distance, the average wind mass density is po = 9 x 10~ 15 cm~ 3 , as inferred 
from the observed H a EW assuming an homogeneous and isotropic wind (C05). 

1 /2 

The correction for dumpiness would require po f (see, e.g., Owocki & Cohen 
2006). The wind /clump speed is fixed to its value at ro, i.e. v o = 1.5x 10 s cm s _1 . 

Here, we are interested on finding the fo dependence of the AiVn, and the 
fo and R c dependences of the number of clumps in the line of sight (N c ) to 
derive, using ANr, restrictions on the wind parameters. N c is relevant given 
the fact that, although a single clump could be optically thick, the probability 
to catch a clump may be significantly smaller than 1. In that case, since the 
XMM-Newton exposure was similar or smaller to the time needed by a clump to 
cover the interclump typical distance ~ R c / fo, we can assume that the duration 
of these observations is not relevant. 

For the simplicity of the employed formulas, we consider the system/X-ray 
emitter /observer geometry for OBS1 and OBS2 to be the same as the one for 
SUPC and INFC. Given the moderate eccentricity of the system, this simplifica- 
tion should not affect much the final results. We also neglect Nr around INFC 
when compared to the value at SUPC, i.e. ANr ~ supc- 

From the model and observational constraints presented above, we can ob- 
tain the formulae for N c and ANr (around SUPC). For N c , we derive first: 



where <j) is the orbital angle from periastron, which is small in our case, and f3 
the angle between the line of sight and the line joining the compact object and 
the star, and (5 = 30° and 60°, which correspond to inclination angles of i = 60° 
(neutron star case) and 30° (black hole case; see C05). From A^h, we can obtain 
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where e is the system eccentricity e = 0.35 (C05), and then 
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Figure 1. The filling factor of the wind as a function of the clump size, 
derived assuming an X-ray emitter located close to the compact object and 
i = 45°, is shown. We note that for reasonable values of i?, c , fo is about 0.01, 
which would imply a mass-loss rate 10 times smaller than that inferred for an 
homogeneous wind. In case R c < 109 cm, /o would not depend anymore on 
Rc, and should be < 0.001. 



The free parameters of our model are the inclination angle i and the wind filling 
factor /o. 

We note that observations impose ANn < 10 21 cm~ 2 . This restriction 
plus the formulae shown above actually imply that, either N c < 1, and then 
fo £ (0.01-0.02) R 2 Ji , or N c > i ) an d then f £ 0.0004-0.003. For N c < 1 and 
a reasonable value of R c , e.g. R c = 0.01 R* = 6.5 x 10 9 cm (R* ~ 6.5 x 10 11 cm; 
C05), the wind filling factor should be fo ^ 0.01. In Fig. [TJ we plot fo as a 
function of the R c for i = 45° under N c < 1, the most stringent condition as 
long as R c > 10 9 . We also plot in Fig. [5] the i-dependence of fo for R c = 0.01 i?* 
(N c < 1). 

3. Discussion 

There are several caveats concerning the adopted wind model. The first one 
is the fact that for i ^ 60° the light rays during SUPC should have a quite 
small impact parameter with respect to the star. It would imply that the wind 
structure should be smoother, the density higher and the velocity lower. This 
would increase AiVn leading to stronger constraints on the real stellar mass- 
loss rate than those derived above. Also, even if we are not in a so extreme 
case as that of % = 60°, the progressive wind acceleration means a lower wind 
speed at distances from the star shorter than rrj, again effectively increasing 
the AJVh- Finally, the clump could have time to expand significantly from 
the region it was formed. Accounting for a clump sound speed ~ 10 6 cm s _1 
(T ~ 10 4 K) and an initial size R c q < 10 10 cm, its size would grow significantly, 
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Figure 2. The filling factor as a function of the inclination angle for a clump 
size of 0.01 i?*. The inclination dependence is rather smooth, and /o becomes 
smaller when i grows. We note that for i > 60° the X-ray would pass very 
close from the stellar surface around periastron/SUPC, implying that /o is 
not meaningful anymore. Actually, for such a i- values, the hydrogen column 
density would be much larger than the one assumed here, as discussed in the 
text. In such a case, soft X-rays should be absorbed even for small M w , and 
they are not. 



i.e. AR C > R c o, increasing the interaction probability and eventually entering 
in the regime described by Eq. ([3]) with the original /n, since it depends only on 

rl/2 

It is worth noting that a wind suffering strong ionization would yield a 
lower AAh than the real value, provided that there are no atoms to be ionized 
by X-rays. It is certainly happening in bright X-ray sources like Cygnus X-l, 
but LS 5039 is a relatively faint X-ray source. The X-ray luminosity around 
SUPC is few times 10 33 erg s , yielding a wind ionization region much smaller 
than the binary system size (B07). 

It could also happen that we are underestimating the real wind absorption 
with our X-ray fits if there were an additional soft X-ray component on the top 
of the power-law one. This possibility cannot be discarded, and may indicate 
the presence of a so far undetected accretion disk. Otherwise, the soft X-ray 
spectrum may be harder, and then we would be overestimating the A Ah- In 
such a case, the constraints on the wind dumpiness and mass-loss rate would 
be even stronger. 

We remark that one can easily derive constraints for the wind properties 
when the emitter is not located right on the compact object but at different 
locations in relation to the star. For that, the angle dependence with [3 is the 
same, but one must modify accordingly the (5 — i relationship. In addition, the 
distance from the emitter to the star center is not the orbital distance anymore, 
and should be replaced by the proper distance to derive Eqs. 1, 2 and 3. This 
allows one for instance to find at which distance the X-ray emitter should be 
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from the star in case the wind were (approximately) homogeneous. Interestingly, 
in such 8i CctSG, ELS already noted by B07, the X-ray emitter should be located at a 
distance larger than the orbital radius from the star and also from the compact 
object. A more detailed study will be presented elsewhere. 

3.1. Consequences of a deep emitter 

Putting aside the possibility that the X-ray emitter is far from the compact 
object, i.e. ^ a from it, there remains only a very porous and therefore weak 
wind as an alternative. The value of the mass-loss rate from the /o values 
derived above is M w ~ 5 x 10 -8 M /yr. This mass-loss rate implies a very 
weak wind momenta and, if the compact object in the system were a non- 
accreting pulsar (see Martocchia et al. 2005; Dubus 2006; Sierpowska-Bartosik 
&; Torres ), the spin-down luminosity of the pulsar wind required to explain 
the observations, of several 10 36 erg s _1 , would generate a shock with a fairly 
wide opening angle, yielding an uncollimated structure (Bogovalov et al. 2008, 
Romero et al. 2007)o This would also mean that the radio pulses would not 
be free-free absorbed during a substantial fraction of the orbit. All this would 
imply that, in the pulsar scenario, the shocked pulsar outflow could hardly be 
the origin of the collimated radio emission, and pulses could be detectable. On 
the other hand, we also note that such a weak stellar wind (/o ~ 0.01) would be 
hardly enough to power the observed non-thermal emission via accretion. We do 
not investigate here neither how a very porous stellar wind can interact with the 
pulsar wind leading to non-thermal emission, nor how a very clumpy wind can 
lead to accretion and jet formation, but these issues should be also addressed. 

Given the mentioned problems in the two proposed scenarios, a far X-ray 
emitter would be preferable, which is not compatible with an X-ray emitter 
between the star and the pulsar, and would favor an emitter based on some 
sort of jet-like structure, possibly like the radio and the TeV emission (see the 
discussion in Bosch-Ramon 2009; Khangulyan et al. 2008, Bosch-Ramon et al. 
2008), regardless the nature of the compact object, as noted in B07. 

4. Summary 

In this work, we show that the lack of significant features due to stellar wind 
absorption in the X-ray spectrum of LS 5039, assuming an emitter close to the 
compact object, indicates that the wind has a filling factor ~ 0.01 or smaller. 
Such a low wind filling factor would imply a stellar mass-loss rate one order of 
magnitude below the values inferred by C05, which would have several implica- 
tions in both the pulsar and the accretion scenario. In the former, such a weak 
wind would imply a very large opening angle for the shocked wind structures. 
In such a case, the lack of radio pulses could not be explained by free-free ab- 
sorption in the stellar wind for a significant part of the orbit, nor the shocked 
pulsar wind would appear collimated unless i were close to 0° or 90°. In the 



2 We note however that opening angles close to 7r would imply an X-ray emitter of a size similar 
to that of the binary system, making our calculations not applicable. 
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accretion scenario, such a low mass-loss rate could hardly sustain the power of 
the observed persistent non-thermal bradband emission. 

An emitter located far from the compact object would not have the prob- 
lems mentioned above for the two proposed scenarios, and the emitter would be 
consistent with a jet-like emitter. 

Acknowledgments. W.B-R. thanks Anabella T. Araudo, Dmitry Khangulyan, 
Stan Owocki and Gustavo E. Romero for fruitful discussion on the topic of this 
work. V.B-R. gratefully acknowledges support from the Alexander von Hum- 
boldt Foundation. V.B-R. acknowledges support by DGI of MEC under grant 
AYA2007-68034-C03-01, as well as partial support by the European Regional 
Development Fund (ERDF/FEDER). 

References 

Aharonian et al. 2005, Science, 309, 346 

Aragona, C, McSwain, V., Grundstrom, E., et al. 2009, ApJ, submitted [astro- 
ph/0902.4015] 

Bogovalov, S. V., Khangulyan, D., Koldoba, A. V., Ustyugova, G. V., Aharonian, F. A. 

2008, MNRAS, 387, 63 
Bosch-Ramon, V., Motch, G, Ribo, M., et al. 2007, A&A, 473, 545 (B07) 
Bosch-Ramon, V., Khangulyan, D., Aharonian, F. A. 2008, A&A, 489, L21 
Bosch-Ramon, V. 2009, A&A, 493, 829 

Casares, J., Ribo, M., Ribas, I., et al. 2005, A&A, 364, 899 (C05) 
Dubus 2006, A&A, 456, 801 

Khangulyan, D., Aharonian, F. A., Bosch-Ramon, V. 2008, MNRAS, 383, 467 
Marti, J., Paredes, J. M., Ribo, M. 1998, A&A, 338, L71 
Martocchia, A., Motch, G, Negueruela, I. 2005, A&A, 430, 245 
Owocki & Cohen 2006, ApJ, 648, 565 

Paredes, J. M., Marti, J., Ribo, M., Massi, M. 2000, Science, 288, 2340 
Paredes, J. M., Ribo, M., Ros, E., Marti, J., Massi, M. 2002, A&A, 393, L99 
Ribo, M., Paredes, J. M., Moldon, J., Marti, J., Massi, M. 2008, A&A, 481, 17 
Romero, G. E., Okazaki, A. T., Orellana, M., Owocki, S. P. 2007, A&A, 474, 15 
Takahashi, T., Kishishita, T., Uchiyama, Y., et al. 2008, ApJ, submitted [astro- 
ph/0812.3358] 

Sierpowska-Bartosik, A. & Torres, D. F. 2008, Astropart. Phys., 30, 239 



